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SECTION  I 


INTRODUCTION 


The  impetus  for  the  present  study  was  a  series  of  roughness  ex¬ 
periments  conducted  in  the  AFWAL  Mach  6  facility.  These  experiments  were 
to  examine  roughness  induced  shock  vorticity  effects  on  boundary  layer 
structure.  Two  types  of  roughness  patterns  will  be  investigated:  a 
square  shape  and  a  diamond  shape.  These  patterns  were  obtained  by  machining 
grooves  in  a  smooth  flat  plate,  as  depicted  in  Figure  1.  The  square  elements 
will  produce  a  detached  shock  resulting  in  variable  vorticity,  while  the 
diamond  elements  will  have  attached  shocks  and  constant  vorticity.  As  of 
the  writing  of  this  report,  only  the  square  shaped  pattern  has  been  tested. 
The  roughness  elements  are  0.02  inches  high  by  0.04  inches  on  each  side 
with  center-to-center  spacings  of  0.08  inches. 


■MILLED  GROOVE 

(NOT  TO  SCALE)  FLOW 


Figure  1.  Square  and  Diamond  Shaped  Roughness  Patterns 


Due  to  the  thin  boundary  layer  at  high  Mach  numbers,  it  is  ex¬ 
tremely  difficult  to  make  detailed  flow  measurements.  Furthermore,  in 
hypersonic  flow  at  high  Reynolds  number,  most  of  the  boundary  layer  is 
supersonic.  In  fact,  even  small  roughness  elements  can  generate  shocklets. 
If  these  shocklets  are  curved,  due  to  the  shape  of  the  elements,  vorticity 
will  be  added  to  the  flow  and  turbulent  quantities  such  as  the  boundary 
layer  growth  rate  will  be  affected.  Because  of  the  need  to  simulate  Mach 
number  and  shock  generated  vorticity  effects  for  rough  surfaces,  one  does 
not  have  the  luxury  of  scaling  up  roughness  elements  as  can  be  done,  for 
example,  in  an  environmental  wind  tunnel. 

On  the  other  hand,  detailed  flow  calculations  can  be  performed 
over  a  rough  surface,  even  for  small  elements  in  thin  boundary  layers.  Such 
parameters  as  Mach  number,  wall  temperature,  and  pressure  gradient  can  be 
easily  simulated.  Profiles,  skin-friction  and  heat  transfer  coefficients, 
and  fluctuating  turbulent  quantities  can  be  calculated.  One  can  also 
examine  different  roughness  element  geometries  with  attached  and  detached 
shocks  to  see  how  vorticity  and  momentum  defects  influence  boundary  layer 
structure. 


Most  roughness  analyses  are  based  on  Nikuradse's  sand-grain  ex¬ 
periments  and  the  "law  of  the  wall"  velocity  profiles  fit  to  the  data. 
Several  correlations  have  been  proposed  to  relate  real  roughness  heights, 
spacings,  and  geometries  to  an  equivalent  sand-grain  roughness  height  so 
Nikuradse's  data  can  be  used.  Examples  of  such  correlations  can  be  found 
in  Betterman  1  ,  Dirling2,  and  White  and  Grabow3.  However,  the  general 
applicability  of  such  correlations  is  unknown.  More  recently,  surface 
roughness  calculations  have  been  performed  by  differential  methods.  Cebeci 
and  Chang"  numerically  solved  the  incompressible  boundary  layer  equations 
employing  an  algebraic  eddy-viscosity  formulation  modified  for  surface 
roughness.  This  modification  was  based  on  Rotta’s  model5,  which  displaces 
the  normal  coordinate  of  the  rough  wall  velocity  profile.  An  expression 
for  this  displacement,  and  the  resulting  mixing  length,  is  given  by  Cebeci 
and  Chang  as  a  function  of  an  equivalent  sand-grain  roughness  height. 


Emphasizing  compressible  flows  for  a  variety  of  edge  and  wall  conditions, 
Hodge  and  Adams6  numerically  solved  the  boundary  layer  equations  and  an 
integral  form  of  the  kinetic-energy-of -turbulence  equation.  Roughness 
effects  were  accounted  for  by  a  form-drag  term  in  the  momentum  equation 
and  by  modifications,  based  on  the  results  of  Healzer  et  al.7,  for  several 
of  the  nine  empirical  constants.  It  is  the  intent  of  this  study  to  utilize 
a  numerical  technique  for  calculating  skin  friction  and  heat  transfer  over 
surfaces  with  real  roughness,  including  the  effect  of  shock  waves,  using  as 
little  empiricism  as  possible. 

Contriving  an  equivalent  sand-grain  roughness  height  for  real 
roughness  heights,  spacings,  and  geometries  is  not  very  satisfying.  A 
physically  more  meaningful  method  that  accounts  for  actual  roughness  effects 
is  that  employed  by  Finson  and  Clarke8,  Lin  and  Bywater9,  and  Christoph 
and  Pletcher10.  These  techniques  calculate  the  form-drag  contributions  of 
individual  elements.  Roughness  elements  are  assumed  to  occupy  no  physical 
space.  The  governing  boundary  layer  equations  are  cast  in  a  form  to  account 
for  the  blockage  effects  of  the  roughness  elements.  Terms  that  act  in  the 
streamwise  direction  are  multiplied  by  [1  -  D(y)/L]  where  D(y)  is  the 
element  width  at  height  y  and  L  is  the  average  center-to-center  spacing. 
Terms  that  act  in  a  direction  normal  to  the  streamwise  direction  are  multi¬ 
plied  by  [  -  D2(y)/L2].  The  effect  of  roughness  is  described  by  a 
sink  term  in  the  momentum  equation  and  by  a  source  term  in  the  static  en¬ 
thalpy  equation. 

When  the  above  approach  is  adopted,  one  must  examine  carefully  the 
turbulence  model  used.  Existing  rough  wall,  mixing-length  models  are  ex¬ 
pressed  in  terms  of  an  equivalent  sand-grain  roughness  height  and  are 
validated  for  boundary  layer  equations  without  blockage  effects  or  source 
and  sink  terms  (e.g.,  Cebeci  and  Chang4  or  Healzer  et  al.7).  Lin  and 
Bywater9  modify  their  turbulent  kinetic  energy  model  equation  to  include 
blockage  effects.  Finson  and  Clarke8  use  a  second-order  closure  approxi¬ 
mation  for  their  turbulence  model,  describing  the  effect  of  roughness  by 
distributed  source  or  sink  terms  in  the  appropriate  equations.  These 


latter  two  methods  are  still  left  with  modeling  constants  whose  roughness 
effects  are  uncertain.  For  the  present  study,  it  was  decided  to  use  the 
analysis  of  Christoph  and  Pletcher10  modified  to  include  shock  wave  effects. 
In  this  analysis  a  two-layer  algebraic  mixing-length  model  is  employed  that 
explicitly  includes  roughness  height,  frequency,  and  type. 

Besides  giving  an  indication  of  measurement  magnitudes  to  be 
expected,  calculations  can  assist  in  designing  experiments.  Roughness  size 
plus  flow  parameters  can  be  determined  to  best  simulate  flight  conditions. 

It  is  hoped  that  computations  and  experiments  progress  hand-and-hand  to 
evolve  an  experimental  program  that  produces  useful  data  for  accurate  tur¬ 
bulence  modeling  and  a  clearer  understanding  of  the  flow  physics. 


SECTION  II 


EQUATIONS  AND  SOLUTION  PROCEDURES 
1.  GOVERNING  EQUATIONS 

For  a  smooth  two-dimensional,  steady,  compressible  turbulent  flow, 
the  standard  boundary  layer  equations  are 

CONTINUITY 


MOMENTUM 


h  <*«>  +  h (p5)  =  0 
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ENERGY 
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The  coordinate  x  is  measured  along  the  surface  from  the  leading  edge  and 
the  coordinate  y  is  perpendicular  to  x.  The  boundary  conditions  imposed 
are  at  y  =  0,  u(x,0)  *  V  (x,0)  =  0,  H(x,0)  *  Hw(x)  or  (3H/3y)y=o  specified, 
and  as  y  +  ®,  u(x,y)  =  ue(x)  and  H(x,y)  =  He(x) .  In  the  above,  v  =  (pv  +  p'v' )/p. 
The  terms  -pv‘u'  and  pv ' h ‘  represent  the  apparent  turbulent  shear  stress  and 
heat  flux,  respectively,  and  must  be  modeled. 


Equations  (1)  -  (3)  are  now  recast,  in  the  manner  proposed  by 
Finson  and  Clarke,  to  account  for  roughness.  As  discussed  in  the  Introduc¬ 
tion,  terms  that  act  in  the  streamwise  direction  are  multiplied  by  [1-D(y)/L], 
and  terms  that  act  normal  to  the  streamwise  direction  are  multiplied  by 
11-  D2(y)/(  L2)].  A  form-drag  term 

-  |  pu2cDD(y)/L2  (4) 

is  added  to  the  momentum  equation,  and  a  source  term 

|  pu3cDD(y)/L2  (5) 


is  added  to  the  static  enthalpy  equation  such  that  the  total  enthalpy  is 
not  altered.  Also,  the  coordinates  are  transformed  by 


and  the  dependent  variables  are  nondimensionalized  according  to 


I  =  £-  ,  V  =  p  v  ( 
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Utilizing  Equations  (6)  and  (7)  and  the  roughness  modifications,  the  con¬ 
servation  equations  become 
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The  functions  B(y)  and  f(y)  are  defined  by 

B(y)  =  1  -  D2(y)/(L2)  (11] 

f(y)  =  (1  -  D(y)/L]/B(y) 

In  the  above,  the  Boussinesq  assumption  has  been  employed  to  evaluate  the 
shear  stress;  that  is,  it  has  been  assumed  that 


3u  — > — r  ,  »  3u  -  3u 

X  =  p  ^  -  pv’u’  =  (M  +  Mt)  ^  =  M 


7. 


It  has  been  further  assumed  that  the  turbulent  heat  flux  can  likewise  be 
represented  through  the  Boussinesq  assumption.  Representing  the  tur¬ 
bulent  conductivity  by 


kt  =  cp  *VPrt  ^13 

where  the  turbulent  Prandtl  number  was  set  equal  to  a  constant  value  of 
0.9,  the  turbulent  heat  flux  can  be  evaluated  as 


I  -pv ' h * 
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(14! 


Equations  (12)  -  (14)  were  employed  to  eliminate  the  Reynolds  stress  and 
heat  flux  terms  prior  to  transforming  the  equations.  The  boundary  condi¬ 
tions  for  the  new  variables  are  given  by 


F(x,0)  =  0,  V(x,0)  =  0 


I(x.O)  =  Iw(x)  or 


specified 


and  as 


D  -*■  »,  F  =  I  *  1.0 


(15] 


In  order  to  include  boundary  layer  edge  changes  for  variable  entropy,  the 
pressure  gradient  is  written  as 


dP  due  .  T  dse 

•  37  ‘  peue  XT  +  peTe  Tx 


(16] 


where  S  is  the  entropy. 


2. 


TURBULENCE  MODELING 


Prandtl's  mixing-length  hypothesis  is  used  to  evaluate  the  apparent 
turbulent  viscosity  according  to 


(17) 


The  mixing  length  is  evaluated  as 

£.  =  0.41  Day  (18) 

in  the  inner  region  and 

s  0.089  6 

in  the  outer  region  where  Da  is  a  modified  van  Driest  damping  function  and 
6  is  the  boundary  layer  thickness.  The  damping  function  and  the  form  of 
the  mixing  length.  Equation  (18),  must  be  specified  for  such  complicating 
effects  as  wall  transpiration  and  surface  roughness.  For  incompressible 
flow  over  an  impermeable  smooth  wall,  van  Driest11  recommended 

Da  *  1  -  exp(-y+/A+)  (20) 

where  A+  =  26.  Pletcher12  modified  this  form  for  transpired  flows  by 
setting 

Da  ■  1  -  exp(-z)  (21) 

where 

2  =  26  (x/iw)1/2(iFT/xw)1/2  (22) 


9. 


and  xpj  is  the  shear  stress  near  the  beginning  of  the  fully  turbulent 
region.  This  form  for  z  was  proposed  to  account  for  shear  stress  varia¬ 
tions  near  a  wall  such  as  occurs  with  wall  blowing  or  roughness. 


If  the  shear  stresses  in  Equation  (22)  are  rough  wall  values,  one 
would  expect  the  smooth  wall  constant,  A+  =  26,  to  change.  Instead  of  try¬ 
ing  to  find  a  new  empirical  constant/function  for  rough  walls,  let  us  try 
multiplying  the  right-hand  side  of  Equation  (22)  by  Twr/iws  where  the  sub¬ 
scripts  r  and  s  represent  rough  and  smooth  values,  respectively.  Then 

z  =  26  (twr/Tws)  (23) 


An  expression  for  tpjr  is  obtained  by  neglecting  the  convective  accelera¬ 
tion  and  by  evaluating  the  momentum  equation  at  the  top  of  the  roughness 
elements  (y  =  K).  Then 


xFTr  Tws  +  2 


i 


K  f (y)  pu  cDD(y) 
B(y)  L2 


dy 


(24) 


where  the  pressure  gradient  term  is  neglected.  As  the  roughness  height 
increases,  the  damping  factor  Da  approaches  one,  as  expected. 

It  is  also  necessary  to  modify  the  mixing-length  formulation. 
Equation  (18),  for  roughness  effects.  Healzer  et  al.7  and  Cebeci  and  Chang4 
suggest  mixing-length  modifications  for  sand-grain  roughness.  For  the 
present  study,  a  mixing-length  was  derived  based  on  the  "law  of  the  wall" 
analysis  of  White  and  Christoph13.  For  illustrative  purposes,  the  de¬ 
rivation  presented  here  considers  incompressible  flow.  First,  approximate 
the  turbulent  shear  using  Prandtl's  mixing-length  hypothesis  according  to 
Equations  (17)  and  (18).  Then,  note  that  near  the  wall  there  is  negligible 
convection  acceleration,  so  that  an  expression  for  tw,  such  as  Equation  (24) 
holds.  Combining  these  results  gives  an  expression  for  Du+/9y+.  For  fully 
turbulent  flow  with  roughness,  one  has 
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3y+  0.41  y+ 

Note  that  for  a  smooth  wall  one  would  obtain 

3u+  _  1 

3y+  0.41  y+ 

where  i  =  0.41  y.  This  suggests  that  for  a  rough  surface 


i  = 


0.41  y 


T  /X 

ws7  wr 


c 


+  vvcDf*y+*  U+  D(y+) 


1/2 


u*  L*  B(y+) 


dy 


(26) 


•  (27) 


For  the  present  study,  the  rough  wall  mixing-length  was  derived  to  include 
compressibility.  The  "law  of  the  wall"  profiles  outlined  above  are  in 
reasonable  agreement  with  the  roughness  channel-flow  data  of  Schl ichting14 . 

These  are  given  in  Reference  10.  This  is  a  partial  verification  for  the 
form  of  the  damping  function  given  In  Equation  (23)  and  the  mixing-length 
modification  just  discussed. 

The  switch  point  from  the  inner  [Equation  (18)]  to  the  outer 
[Equation  (19)]  model  was  made  In  accordance  with  the  low  Reynolds  number 
modification  presented  in  Reference  15.  This  provides  that  the  switch  be 
made  at  the  point  where  first  becomes  equal  to  t0  provided  that  y+  >.  50. 

If  y+  <50  when  t-j  s  Iq,  then  the  switch  point  is  delayed  until  y+  becomes 
equal  to  50  and  i0  becomes  equal  to  £]  at  y+  =  50.  This  prevents  the 
suppression  of  the  fully  turbulent  portion  of  the  velocity  profile  that 
can  occur  at  low  values  of  ReQ.  This  effect  persists  to  higher  and  higher  vaK.es 
of  Reft  as  the  Mach  number  increases. 


3. 


TREATMENT  OF  SHOCKS 


In  a  supersonic/hypersonic  boundary  layer  it  is  possible  to 
generate  shocklets  even  for  small  roughness  heights.  The  analysis  of 
Christoph  and  PI  etcher 10  was  modified  to  account  for  the  effect  of  shock- 
lets.  To  the  author's  knowledge,  the  only  data  available  showing  shocklets 
were  obtained  by  A.  Flore  of  AFWAL  .  These  data  for  the  square  pattern  of 
Figure  1,  form  the  basis  of  the  present  modifications.  Photographs  of  the 
shocklets  are  not  shown  in  this  report  because  they  would  not  reproduce 
clearly.  Qualitatively,  the  shocks  start  at  the  sonic  line  and  are  detached, 
and  they  approach  the  flow  Mach  angle  at  the  boundary  layer  edge.  The  axial 
extent  of  the  shock  before  merging  with  the  boundary  layer  edge  is  approxi¬ 
mately  equal  to  the  roughness  element  spacing.  A  reproduced  drawing  of  the 
Schlieren  is  shown  below.  It  is  impossible  to  describe  the  shock  shapes 
accurately  from  the  available  photographs.  Instead,  the  shock  shape  is 
calculated  by  assigning  values  of  M  sin  0  (code  input  parameter)  where  0 
is  the  shock  angle. 
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Figure  2.  Shocklets  off  Roughness  Elements 


It  is  believed  that  Msin9  is  slightly  greater  than  1.0,  i.«  .,  shocks  art 
weak.  The  reason  for  this  belief  is  that  calculations  made  without  shocks 
are  in  reasonably  good  agreement  with  profile  data  (discussed  in  Results 
section  of  this  report).  Also,  from  the  photographs  the  shock  shape  appears 
invariant  in  the  flow  direction.  In  this  study,  the  shock  shape  is  kept 
similar  at  all  axial  stations,  expanding  with  the  growing  boundary  layer. 

For  determination  of  the  initial  shock  shape,  no  shocks  are  allowed  over 
the  first  few  elements  so  that  a  rough  wall  Mach  number  distribution  through 
the  boundary  layer  can  be  determined.  Finally,  oblique  shock  relations  are 
used  to  calculate  flow  variables  behind  the  shocklets.  Because  of  the  un¬ 
certainty  in  the  shock  angles  at  the  edge  of  the  boundary  layer,  edge 
variations  of  velocity  and  entropy  have  been  set  to  zero  for  the  results 
presented  in  this  report. 

4.  METHOD  OF  SOLUTION 

Equations  (8)  -  (10)  are  solved  numerically  by  a  fully  implicit 
finite-difference  procedure.  As  shown  in  Figure  3,  non-uniform  grid  spacing 
in  both  the  x  and  y  directions  is  used. 
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Values  are  known  at  station  i  and  unknown  at  downstream  station  i  +  1.  Tte 
equations  are  solved  in  an  uncoupled  manner  in  the  sequence:  momentum,  con¬ 
tinuity,  and  energy.  The  finite-difference  formulation  used  depends  on 
the  x  station.  At  station  1  there  is  no  x  variation.  Central  differences 
are  used  for  the  y  variations  at  all  stations.  For  example. 


3y  (Ay++  AyJ1 


«#:  -  F]+1> 


Ay. 

Ay+ 


+  (F 


i+1 


-  Fi-i) 


Ay+ 

Ay- 


(28) 


where  V  denotes  the  most  recent  evaluation  of  V  at  i+1  level .  Unequal  grid 
spacing  in  the  normal  direction  is  implemented  by  a  geometric  progression 
such  that  the  ratio  of  two  adjacent  normal  coordinate  increments  is  a  con¬ 
stant.  This  ratio  was  set  at  1.08  for  the  calculations  reported  here.  For 
the  general  continuing  calculation  we  want  to  use  information  from  the  two 
previous  stations  to  obtain  second-order  accurate  representations  of  the 
x-derivative  terms  at  the  i+1  level  and  to  obtain  a  second-order  accurate 
representation  of  the  velocities  that  appear  as  coefficients  in  the  con¬ 
vective  terms.  However,  for  the  second  calculation  station  we  only  have 
information  at  one  station  so  we  use  a  more  conventional  fully  implicit 
scheme.  At  x  station  2,  for  example. 


xF  — 
*  dx 
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Ax, 
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a 

where  Fj  values  are  taken  as  known  values  at  the  i  level  initially  and 
updated  if  iterations  are  used.  It  is  generally  not  crucial  that  they  be 
updated.  For  all  other  calculation  stations  beyond  the  second,  streamwise 
derivatives  are  approx  lira  ted  to  second-order  accuracy  by  using  three-point 
difference  representations.  Second-order  accuracy  for  the  convective  terms 
is  maintained  without  Iterations  by  the  use  of  extrapolated  values  of  the 
coefficients.  For  example. 


xF  =  x1+1  UxT  (Dx3F.1+1  -  Dx2F -1  +  DxlF.1-1  ) 

a  X  J  J  J 


(30) 
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Dxl  =  Ax+/(Ax_  (Ax+  +  Ax_) ) 

Dx2  =  (Ax+  +  Ax_)/(Ax+  +  Ax_) 

Dx3  =  (Ax_  +  2Ax+)/(Ax+(Ax+  +  Ax_) ) 


and 
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Note  that  is  the  above  difference  representations,  F2  has  been  linearized 
by 


F2 


i+1 


(33) 


Each  difference  equation  for  the  momentum  and  energy  equations 
has  three  unknowns  (Fj*j,  f1+1,  Fj*j)  ancl  can  be  arranged  in  the  form 


A  F^+^  +  B  F^+^  +  C  F^+^  =  D. 
\j-l  Yj  Uj 


(34) 


This  system  of  equations  is  a  tri -diagonal  form  and  is  solved 
efficiently  by  the  Thomas  algorithm.  It  is  important  to  keep  Bj >  0  and 
Aj,  Cj  <  0  or  unphysical  results  with  wiggles  can  follow.  If  Aj >  0,  the 
central  difference  representation,  e.g..  Equation  (28),  is  replaced  by  a 
forward  difference 
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and  if  C . >  0,  a  backward  difference 
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is  used.  The  continuity  equation  can  be  solved  explicitly  for  Vi 
Several  features  of  the  present  difference  procedure  are  as  described  by 
Harris1*.  For  the  final  closure  of  equations,  the  density  is  obtained  from 
the  ideal  gas  relations 


pv/p  s  T/Tw’  h  =  cpT 

and  the  molecular  viscosity  from  Sutherland's  formula 

1/,?  lb  sec 
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u  =  2.27  x  10  8(x  +  {98>6/T) 
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(38) 


The  skin  friction  coefficient  is  calculated  from  the  smooth  wall 
contribution  plus  the  form-drag  contribution 
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Heat  transfer  coefficients  cannot  be  calculated  in  a  similar  manner  be¬ 
cause  there  is  no  heat  transfer  mechanism  analogous  to  form  drag.  A 
simple  technique  suggested  by  Finson  and  Clarke  was  adopted  in  this  study. 
That  is,  the  heat  transfer  augmentation  due  to  roughness  is  equal  to  the 
square  root  of  the  skin  friction  augmentation  due  to  roughness.  Physically, 
this  is  reasonable  because  velocity  fluctuations  are  increased  by  rough¬ 
ness  but  temperature  fluctuations  are  hardly  changed  by  roughness  and 
tw  -  u  V  and  qw~  vT". 
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SECTION  III 


DISCUSSION  OF  RESULTS 

Comparisons  were  made  to  the  smooth  and  square  roughness  pattern 
data  taken  by  Fiore  in  the  AFWAL  Mach  6  facility.  Profile  data  were  ob¬ 
tained  at  17.15  inches  from  the  test  plate  leading  edge  for  tunnel  total 
pressures  of  approximately  700,  1400,  and  2000  psia  and  a  tunnel  total 
temperature  of  1100°R.  Both  probe  and  laser  velocimeter  measurements  were 
taken.  Surface  pressures  and  temperatures  were  also  measured  along  the 
plate  length.  A  summary  of  the  rms  edge  and  wall  conditions  for  the  smooth 
and  rough  wall  runs  are  given  in  Table  1. 


TABLE  1,  NOMINAL  TEST  CONDITIONS 


CONFIG. 

X 

IN. 

\ 

Po 

psia 

°ToR 

Pw 

psia 

k 

Pe/K, 

ft’1 

Rex 

SFP 

17.20 

5.92 

1979.78 

1086.54 

1.36 

601.90 

2.521 x  107 

3.613 x  107 

SFP 

5.89 

1402.54 

1078.39 

1.00 

635.51 

1 .829  x  107 

2.621 x  107 

SFP 

5.89 

703.29 

1067.65 

0.50 

648.69 

0.9309 x  107 

1.334 x  107 

RFP 

5.77 

2001.03 

1131.85 

1.61 

691.94 

2.557 x  107 

3.665  x  107 

RFP 

5.76 

1799.68 

1125.33 

1.47 

751.98 

2.329  x  107 

3.338 x  107 

RFP 

\ 

5.74 

1402.99 

1120.76 

1.17 

778.89 

1.843 x  107 

2.642 x  107 

RFP 

5.67 

704.27 

1114.24 

0.63 

789.64 

1.008 x  107 

1.445 x  107 

(SFP  =  SMOOTH  FLAT  PLATE,  RFP  -  ROUGH  FLAT  PLATE) 


Calculations  were  performed  for  all  three  tunnel  pressures.  Since 
results  are  similar,  profile  comparisons  are  shown  only  for  P0  =  1400  psia. 
Figure  4  is  a  plot  of  the  smooth  and  rough  wall  velocity  profiles  obtained 
with  standard  probes  and  a  laser  velocimeter  (L.V.)  method.  Both  experi¬ 
mental  techniques  produce  profiles  that  agree  very  well  with  each  other. 

This  agreement  applies  to  both  the  smooth  and  rough  models.  It  should  be 
noted  that  the  velocity  profile  for  the  rough  flat  plate  is  not  as  full  as 
that  for  the  smooth  flat  plate.  This  is  caused  by  a  momentum  loss  in  the 
flow  field  around  the  roughness  protuberance  resulting  from  both  local  shock 
waves  and  wakes  generated  by  the  protuberances  themselves.  Predictions  are 
in  good  agreement  with  the  data  for  both  the  smooth  and  rough  surfaces. 

Symbols  represent  the  data  and  curves  the  theory.  A  value  of  Cp  =  0.6  was 
used  for  all  rough-wall  calculations  presented  in  this  report. 

Static  temperature  profile  predictions  are  shown  in  Figure  5.  For 
a  constant  value  of  y/6  the  rough  surface  static  temperature  ratio  is  greater 
than  the  smooth  wall  case.  For  example  at  y/6  =0.275  the  ratio  is  3.45  for 
the  rough  wall  as  compared  to  2.50  for  the  smooth  wall  case.  The  smooth  wall 
computations  agree  reasonably  well  with  the  static  temperature  ratio  measure¬ 
ments  for  0.3  s  y/6<  1.0;  however  in  the  range  for  0  <  y/6  <  0.3  the  compu¬ 
tations  underpredict  the  measured  values  of  T/Te.  This  may  be  a  problem 
with  locating  the  temperature  edge  value  in  the  computer  code.  Temperature 
results  near  the  wall  are  sensitive  to  the  procedure  for  defining  the  tempera¬ 
ture  edge  conditions.  Rough  wall  results  look  good  but  this  may  be  fortuitous 
in  light  of  the  smooth  wall  results.  It  should  be  noted  that  very  near  the 
wall  the  theory  indicates  a  peak  in  T/Te.  This  peak  was  not  noted  in  the 
experimental  data  simply  because  the  probe  size  did  not  permit  obtaining 
data  in  the  region  0  <  y/6  <  0.07.  The  peak  predicted  for  the  rough  wall 
may  be  too  high,  because  an  increase  in  the  rough  wall  mixing  length  near 
the  wall  for  the  momentum  equation  results  in  a  corresponding  increase  in 
the  temperature  calculated  from  the  energy  equation.  However,  for  a  rough 
surface  heat  transfer  around  the  roughness  elements  is  principally  by 
molecular  conduction. 


Underprediction  of  the  smooth  wall  static  temperature  near  the 
plate  surface,  along  with  slight  discrepancies  in  the  velocity  predictions, 
results  in  the  total  temperature  being  underpredicted.  This  is  shown  in 
Figure  6.  Again,  rough  wall  values  look  good.  Basically  the  local  total 
temperature  measurements  are  the  same  for  both  the  smooth  and  rough  flat 
plates.  The  computations  also  indicate  that  the  temperature  distribution 
is  relatively  invariant  with  roughness.  The  local  Mach  number  distribution 
is  presented  in  Figure  7.  In  general  the  local  Mach  number  for  the  rough 
surface  is  lower  than  the  smooth  wall  case.  This  is  believed  to  be  due  to 
a  decrease  in  momentum  caused  by  the  roughness  protuberances.  For  both  the 
smooth  and  rough  walls,  the  Mach  number  in  the  outer  part  of  the  boundary 
layer  is  underpredicted  because  the  static  temperature  is  overpredicted  and 
the  velocity  is  slightly  underpredicted . 

The  Pitot  pressure  through  the  boundary  layer  is  presented  in 
Figure  8.  The  measurements  indicate  that  the  local  Pitot  pressures  for 
the  rough  flat  plate  are  everywhere  less  than  the  smooth  flat  plate  case. 

The  computations  agree  reasonably  well  with  the  measurements  for  both  the 
smooth  and  rough  flat  plates,  except  in  the  outer  region  where  the  Mach 
number  calculations  disagree. 

Figure  9  is  interesting,  but  not  as  dramatic  as  results  for 
P0  =  2000  psia.  The  smooth  and  rough  wall  data  collapse  on  the  same  curve 
when  T/Te  is  plotted  against  u/ue.  (Note  that  the  differences  shown  in 
Figure  9  result  from  very  different  values  of  Tw/Te  for  the  smooth  and 
rough  wall  runs.)  This  implies  that  a  Crocco  type  approximation  is  valid 
for  rough  surfaces,  even  at  Mach  6. 

Skin  friction  coefficients  were  obtained  by  Fiore  using  both  a 
Preseton  tube  and  by  iterating  on  "law  of  the  wall"  velocity  profiles.  In  the 
smooth  wall  case  the  Preston  tube  values  were  used  while  the  iterative  skin 
friction  coefficients  for  the  rough  wall  are  reported.  These  results,  along 
with  predictions,  are  shown  in  Table  2.  Predictions  are  in  very  good  agreement 
with  the  data,  except  for  the  rough  surface  at  P0  =  2000  psia.  The  data  show 
an  increase  in  c*  in  going  from  P0 =  1400  psia  to  P0 =  2000  psia  while  the 
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TABLE  2.  SKIN  FRICTION  COEFFICIENT  COMPARISONS  AT 
x  =  17.15  INCHES 


CONFIGURATION 

Me 

Po 

psia 

k 

Jw 

°R 

cf  x  103 
(Experiment) 

Cf  x  103 
(Theory) 

SFP 

5.72 

701.45 

1080.03 

604.62 

1.202 

1.131 

SFP 

5.80 

1398.72 

1088.39 

604.62 

1.002 

1.000 

SFP 

5.94 

1995.38 

1098.60 

664.62 

0.898 

0.910 

RFP 

5.44 

705.61 

1098.49 

764.37 

1.654 

1.543 

RFP 

5.60 

1405.19 

1106.74 

754.18 

1.347 

1.368 

RFP 

5.53 

2003.53 

1122.94 

654.52 

1.622 

1.363 

predictions  show  a  leveling  off.  It  should  be  noted  that  the  experimental 
rough  wall  Cf  values  may  not  be  very  accurate  because  of  the  dearth  of  data 
in  the  log  portion  of  the  profiles. 

It  should  also  be  mentioned  that  the  predicted  smooth  wall  values 
of  the  boundary  layer  thicknesses  are  in  good  agreement  with  the  experi¬ 
mental  data.  While  both  the  measured  and  calculated  rough  wall  boundary 
layer  thicknesses  were  found  to  be  essentially  independent  of  Reynolds  number, 
the  magnitudes  were  quite  different  (6exp  =  0.38  inches  and  6theory=0‘52  inches) 

Finally,  in  the  calculations  reported  above,  shocks  were  assumed 
weak  (M  sin  8  =  1.025).  This  gave  shock  shapes  that  appeared  similar  to 
Schlieren  obtained  by  Fiore.  Calculations  were  also  made  without  shocks. 

There  was  no  discernible  effect  of  the  shocklets  for  Mach  number,  total 
temperature,  and  total  pressure  distributions.  Near  the  wall,  slight  shock 
effects  were  noticed  for  velocity  and  static  temperature  distributions. 

With  the  shocks,  u/ue  decreased  less  than  4  percent  for  y/6  <0.2  and  T/Te 
increased  less  than  2  percent  for  y/6  <  0.5.  Outside  of  these  y/6  ranges, 
u/ue  and  T/Te  were  essentially  unchanged.  For  the  case  of  tunnel  total  pres¬ 
sure  of  1400  psia,  modeling  the  shocks  in  the  calculations  increased  the 
boundary  layer  thickness  by  2.5  percent  and  decreased  the  skin  friction  coef¬ 
ficient  by  1  percent.  26. 


SECTION  IV 


CONCLUS I ONS/ RECOMMENDAT IONS 

The  fully  implicit  finite-difference  boundary  layer  procedure  of 
Christoph  and  PI  etcher 10  for  flow  over  rough  surfaces  has  been  extended  to 
include  the  effect  of  shocks  off  the  roughness  elements.  The  method  employs 
a  two-layer  algebraic  mixing-length  model  that  explicitly  accounts  for 
roughness  height,  spacing,  and  geometry.  The  equivalent  sand-grain  concept 
is  not  used.  Oblique  shock  relations  are  used  to  calculate  flow  conditions 
through  the  boundary  layer  behind  each  shock. 

Computations  are  presented  for  smooth  and  rough  wall  data  obtained 
at  AFWAL  .  In  general,  the  agreement  between  experiment  and  theory  is  very 
good.  The  smooth  wall  static  temperature  predictions  are  in  disagreement 
with  the  data  near  the  surface.  It  is  felt  that  this  could  be  a  problem 
with  the  code  search  procedure  for  Te,  and  not  with  turbulence  modeling. 
Also,  for  both  the  smooth  and  rough  plates,  the  theory  predicts  a  maximum 
temperature  near  the  wall  that  is  not  captured  by  the  measurements.  Smooth 
wall  skin  friction  coefficients  agree  very  well  with  the  data.  The  rough 
wall  skin  friction  predictions  are  encouraging,  even  though  the  accuracy  of 
the  skin  friction  coefficients  obtained  from  the  rough  wall  profiles  is 
uncertain. 


Several  extensions/improvements  to  the  method  described  in  this 
report  are  suggested.  These  are: 

(1)  Higher  order  turbulence  models  should  be  used.  Possibilities 
are  a  lag  equation  model  and  a  turbulence  kinetic  energy  - 
turbulence  dissipation  rate  (k-e)  model. 

(2)  Mass  addition  should  be  added  as  a  boundary  condition  and 
in  the  turbulence  model  so  that  coupled  roughness/blowing 
effects  can  be  examined. 


(3)  From  other  data  comparisons,  it  has  been  observed  that  the 
code  of  Christoph  and  PI  etcher10  does  not  accurately  pre¬ 
dict  rough  wall  heat  transfer.  The  reason  is  that  an  in¬ 
crease  in  the  mixing-length  for  the  momentum  equation 
results  in  an  increase  in  the  heat  transfer  calculated  from 
the  energy  equation.  However,  for  a  rough  surface  heat 
transfer  near  the  wall  is  principally  by  molecular  conduc¬ 
tion.  A  conduction  sublayer  model  within  the  framework  of 
the  above  analysis  is  needed. 
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